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Abstract 

In this article, we show how to recast the Hubbard model in one 
dimension in a hydrodynamic language and use the path integral ap- 
proach to compute the one-particle Green function. We compare with 
the Bethe ansatz results of Schulz and find exact agreement with the 
formulas for spin and charge velocities and anomalous exponent in weak 
coupling regime. These methods may be naturally generalized to more 
than one dimension by simply promoting wavenumbers to wavevectors. 

1 Introduction 

The Hubbard model is one of the most extensively studied models in Condensed 
Matter Physics. It it the simplest example of an interacting Fermi system on 
a lattice. In one dimension the ground state and excited states may be written 
down explicitly using Bethe ansatz[l]. The collection of reprints by Korepin and 
Essler[3] is particularly useful. The pioneering work of Schulz [2] in computing 
the spin and charge velocities and the anomalous exponent for all values of 
the onsite repulsion U is quite significant. There have been other attempts 
notably by Weng et.al.[4] in recasting the Hubbard model in a path integral 
language. Unfortunately all these methods are restricted to one dimension. One 
notable exception is the work by Liu [5] that uses eigenfunctional theory. Our 
approach is quite different from all these and uses the hydrodnamic formulation 
recently developed by the author. It it simple and powerful as we shall see 
since it is naturally generalisable to more than one dimension and yields the 
important correlation functions quite easily. Our approach is able to give the 
right functional dependence of the various quantities such as spin and charge 
velocities and exponents (as a function of U) for small U. The formalism is 
naturally generalisable to more than one dimension by simply promoting wave 
numbers to wavevectors. Unfortunately we are unable to probe the large U case 
which is relevant to high Tc since the formalism does not reproduce the right 
results in one dimension. 
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2 The Hydro dynamic Formulation 



The program of quantizing hydrodynamics has a long and distinguished his- 
tory. Landau [6] and his students were among the first to attempt this. Later 
on Sunakawa et.al. [7] and others - notably Rajagopal and Grest [8] took this 
program further. Dashen, Sharp, Menikoff and Goldin[9] in the seventies intro- 
duced many of these ideas. Recently, Jackiw and collaborators[10] have revived 
interest in this approach in the context of rclativistic quarks. In our earlier 
work, we introduced the density phase variable ansatz for fermions[ll]. We also 
note that Rajagopal and Grest [8] had already in the seventies pointed out the 
need for having a nonzero-phase functional found in the density phase variable 
ansatz. In our earlier work [11] we made a first pass at computing the phase 
functional. This attempt yielded an answer that in retrospect is quite wrong. 
Upon closer examination the ?7o(q) of our earlier work[ll] is imaginary when it 
was postulated to be real(for small q). So far the author has avoided this issue 
by taking refuge under the the sea-boson approach that enables us to derive 
the momontiim distribution, anomalous exponents, quasiparticlc residue and so 
on without yielding the full dynamical propagator which is of interest only be- 
cause it contains information about quantities just mentioned. If one is able 
to compute them without having to compute the full propagator so much the 
better. However, there are physical problems in which the full propagator is 
important. The X-ray edge problem[14] is one such. In fact we tried using the 
DPVA to compute the X-ray edge spectra in a preprint[13] and found that we 
obtain the right answers in one dimension but the answers in higher dimensions 
were inconsistent with Mahan's exact results[14]. 

In an earlier preprint, after much reflection, we chose to dismiss the ap- 
proach that only uses the hydrodynamic variables namely the density and its 
conjugate as 'myopic' (mypoic bosonization) . This is beacause a hamiltonian 
formulation in terms of the hydrodynamic variables is unable to distinguish be- 
tween fermions and bosons. We have to further decompose these variables in 
terms of linear combination of oscillators in order to distinguish between the two 
statistics. However, the sea-boson approach is not without its share of prob- 
lems. For one it docs not generalise to finite temperatures easily. Also the full 
dynamical propagator is not reducible to quadratures due to a technical difii- 
culty. Both these problems may be resolved in an approach that incorporates 
only the hydrodynamical variables. We show in this preprint, that the path 
integral approach is an avenue to distinguish between the statistics when using 
only the hydrodynamical variables. This approach that only uses hydrodynamic 
variables was not expected to work out for fermions, since one has to take into 
account the extended nature of the Fermi surface. However, this idea seems too 
important to pass up. In particular, the natural manner in wliic;h gauge theory 
may be studied in this approach[15] makes this effort for fermions worthwhile 
and urgent. 

In this section, we introduce the hydrodynamic formulation that has been 
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developed by the author. It involves writing the field variable in trems of ob- 
servables such as currents and densities. In the long-wavelength limit, it can 

be shown that current algebra (that is, the mutual commutation rules between 
currents and densities) is obeyed only if the current operator is expressible as 
shown below. 

j(x,i) = -p(x,i) vn(x,i) (1) 

Here n is the potential for the velocity and is conjuugate to the density p. Thus 
in the hydrodynamic limit the velocity operator is irrotational. The hydrody- 
namic part of the field variable may be written in a polar form. 



^slo^.• (x, t) = eiA(W;x,t)e-ffl(x,t) ^^^f^^ (2) 

For fermions, we expect ipsiow (x, t) to be a Grassmann variable. But we shall 
take the point of view that the fermionic nature of the field is captured at the 
level of the propagator by the introduction of the phase functional A([p];x,t). 
The fermionic KMS boundary conditions is obeyed since a position indepen- 
dent global Klein factor is able to capture this [12]. Thus the program involves 
expanding A{[p]:x.,t) in powers of the density fluctuations and making contact 
with the free theory and fixing the coefiicients. The rest of the discussion is 
similar to our earlier work on bosons [15], therefore we shall not dwell on those 
details. Suffice it to say that if we expand the action in powers of the desnity 
fluctuations and retain only the harmonic parts, in order to recover the right 
desnity-density correlation functions, we have to set, 

H[p]-<in) = C(qn) p_q _„ (3) 

and A{[p];x,t) = X^q_„ e"i''e^"*A([p]; q, n). The coefficient may be computed 
as follows. 

and {pci^nP-ci,-n)o is the density-density correlation function of the free the- 
ory obtained from elementary considerations. In one dimension, for spinless 
fermions we may write, 

where vp is the Fermi velocity and Zn = 27rn//3 is the bosonic Matsubara 
frequency. 
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3 Hubbard Model 



The Hubbard model in one dimension [14] may be written down as follows(here 
G = 27r/a). 

H = -2tJ2 COs{ka) 4^Cfe<,+ ^ ^ Pq^P-ql + ^ ^ Pq^P-q+Gl + -^ X] Pl^P-Q-Gl 

ka ° g " g " q 

(6) 

First we replace the cosine dispersion by a parabolic one. efe = —2t cos{ka) = 
Co + cifc^. From this we find cq = —2t and we require that the slope of the 
dispersion at fc = ifci? be identical to the cosine dispersion. This means c\ = 
{talkp)sin{kpa) = l/(2m). 

Thus we may write a quadratic action for the one band Hubbard model in 
one dimension, including umklapp processes in the hydrodynamic language. 

^Hubb ~ ^ ^ ( '^P^n)Pqa,n-^q(T,n ~t~ ~ ^ ^ ^ 2771 
qa,n qtT,n 

+if3 X Z„C{q7l) pqa,nP-qa-n + ^ PqhnP-ql-n 

aq^On q^On 

if3U \ ^ ipU \ ^ 

+ 770" 2^ Pqt,nP-q+Gl-n + 2^ Pg|,„p_g_Gi,-„ (7) 

q^On q^On 

where G = 2-^1 a. The quadratic action implicitly ignores three body density 
correlation functions. It is not clear why this is valid except that it renders 
the path integrals tractable. However wo may expect to find nontrivial results 
already at the harmonic level with umklapp processes for strong coupling. Since 
we are considering umklapp process which involves large momentum transfer, 
we have to make sure that the C{qn) is evaluated in general for both small and 
large q. In our earlier work we showed, 

- nPq.,nP-q.,-u\ - -^NW- - 

where is the total number of electrons including both spins and (pgcr,nP-gcr,-ra)o 
is the density- density correlation function of the free theory evaluated using el- 
ementary considerations. The field variable is now given by, 

Here tpsiow is the hydrodynamic (slow) part of the field. Also 77 is a 'fudge factor' 
needed to make sure that we recover the right exponents. It is a numerical factor 
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like 1/2 or 2. In the first instance, we may neglect umklapp processes. This 
corresponds to weak coupling and small q. 
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9." 2A„ „JVU'^9-I-''' 



NO 



/(4Ag,„) 

y (^) (/3^n)' + (/3^n) C(g,n) 



2PN' 



C{q,n) = 



4eF 



If we set 7/ = 2 we find, 



(10) 



Q _ \ ( p-iqx „-Znt _ -iqx „-z„t \ 

~ pzn \ e e J 



E„.„ ^ 



q(x-x )„-«„(*-* )_ ) «„( 



(11) 



For weak coupling we have, 

_2tU±Y^ 1 r°° dq / )-Vpj(t-t' )]_ iql{x-J ) + Vpjq{t-t' )]\ 



n 



(12) 



1 

Here, = i>f ^1 ± ^ ^ • The anomalous exponent related to the momen- 



tum distribution is given by, 



7 = 2a = 



16(fc2/TO)2 



(13) 



The Bethe ansatz result of Schulz shows that 7 = U'^ / [Att'^v^.). In units such 
that a = 1 near half filling we have kp = 7r/2 and m = 7r/4. Thus in our case 
we have, 7 = U'^j^. Prom the result of Schulz also we find, 7 = U'^ / {IGTr"^) . 
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Therefore this approach gives the right spin and charge velocities (especially 
close to half- filling), and the right anomalous exponent for small U. 

However, for large U, this approach does not give us the right qualitative 
behaviour as it predicts 7 ~ VU. According to the Bethe ansatz solution the 
anomalous exponent saturates to a value 7 = 1/8. For strong coupling, we are 
unable to make progress. We have tried using the no-double occupany constraint 
since this is easy to implement using the hydrodynamic formulation, but we have 
been unable to make it work. The umklapp terms are also unfortuantely of little 
use. Thus we shall not advocate the use of this approach for strong coupling. 
In two spatial dimensions we expect to find that the system is a Landau Fermi 
liquid. We shall not carry out this calculation since it is not very interesting. The 
main purpose of this article is to highlight the usefuUness of the hydrodynamic 
approach and its generality so that in future publications we may use this to 
study disordred systems and the like. 
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